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ABSTRACT 


Values  of  the  stability  derivatives  C  and  C  have  been 

.  “a  mq 

determined  experimentally  for  cones  with  semi-vertex  angles  of  45°  and 

9°,  and  the  results  compared  with  theoretical  estimates  and  with  other 

experimental  data.  C  is  seen  to  be  a  function  of  angle  of  attack. 

ma 

and  C  is  found  to  be  positive,  Indicating  negative  damping.  The 

“q 

Influence  of  viscous  effects  upon  these  results  is  discussed. 

The  experiment  was  conducted  with  freely  oscillating  models 
in  hypersonic,  low  density  flow,  and  is  one  of  the  first  dynamic  sta¬ 
bility  tests  performed  under  these  flow  conditions.  C  deviates 

ma 

sharply  from  linear  theory,  the  non-linearity  being  even  stronger 
than  predicted  by  durrent  non-linear  Newtonian  flow  theory. 

C  was  found  to  be  positive  in  sign,  meaning  negative 

mq 

damping  and  an  input  of  energy  to  the  model  from  the  flow.  This 
instability  is  a  marked  deviation  from  the  ordinary  assumption  of 
positive  stability  and  certainly  warrants  further  study  and  inves¬ 
tigation.  The  cause  of  the  instability  remains  undetermined;  the 
effects  of  viscous  flow  and  shedding  of  vortices  is  discussed  in 
this  connection. 
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LIST  OF  SYMBOLS 


(Other  Symbols  are  defined  as  used) 


b 


L 


M 

Re/ln 

S 


t 

T 


u,  w 


coefficient  of  damping  term  ■  (-c^  ) 

<1 

standard  non-dimensionalized  moment  coefficient 

restoring  moment  derivative 

damping  derivative 

average  values.  See  Section  5 

rate  of  decay  of  oscillation  (°/sec) 

mass  moment  of  inertia  about  axis  of  oscillation 

2 

coefficient  of  restoring  moment  term  «*  ^  (-C  ) 

ma 

reference  length  ■  axial  length  of  cone 
Mach  number 

Reynolds  number  per  inch  »  ^ 
reference  area  *»  base  area  of  cone 

time 

period  of  oscillation  ■  — 

r  05' 

linear  flow  velocity  components  parallel  to  and  perpendicular 
to  axis  of  revolution,  respectively 


iv  - 


u 


velocity  of  undisturbed  airflow 


X 

a 

a 

7 

P 

<P 

9 

<x> 


distance  from  nose  of  cone  to  center  of  oscillation 
angle  of  attack 

oscillation  angle  or  envelope  of  angle  of  attack 
cone  semi-vertex  or  half  angle 
mass  density  of  atmosphere 
phase  angle 
angle  of  rotation 
angular  velocity  or  frequency 


-  v  - 


1 


1.  INTRODUCTION 

The  motion  of  an  oscillating  cone  is  described  by  the  well- 
known  second  order  differential  equation  containing  inertia,  damping 
and  restoring  moment  terms.  These  terms  can  be  evaluated  experimentally 
by  study  of  the  motion  in  a  wind  tunnel  airstream.  This  report  covers 
such  a  test  in  low  density,  hypersonic  flow;  it  is  in  the  nature  of  a 
feasibility  study  since  there  is  little  background  of  dynamic  stability 
tests  under  such  conditions. 


1.1  Theoretical  Basis 

The  linearized  equation  of  motion  of  natural  short  period 
longitudinal  oscillations  is  given  by  Laitone*  as 
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This  is  an  ordinary  second  order  linear  differential  equation,  with  the 
linearization  implying  small  values  of  0. 


1.2  Induced  Oscillation  Technique 

It  was  first  proposed  (by  Professor  E.  V.  Laitone  of  the 
University  of  California)  that  the  stability  coefficients  of  a  cone  be 
evaluated  by  an  Induced  oscillation  technique  in  the  wind  tunnel.  The 
cone  model  would  be  fixed  on  a  shaft  (sting)  passing  through  and  per¬ 
pendicular  to  the  axis  of  revolution  of  the  cone.  This  sting  would  be 
mounted  in  bearings  and  the  model,  with  sting,  would  be  free  to  revolve 
in  these  bearings.  However,  the  bearing  supports,  and  therefore  the 
model,  would  be  oscillated  in  a  plane  perpendicular  to  the  flow,  with 
a  velocity 

z(t)  «  aco  sin  cot  (3) 

where  a  =  amplitude  of  oscillation 

z  ■  linear  displacement  of  sting 
co  «  frequency  of  oscillation 

Then  in  a  wind  tunnel  of  velocity  U  **  u 


w 

U 


9  +  sin  cot 


and  Equation  (2)  becomes 

0  +  b0  +  k0  ■  -  k  — j  sin  cot 
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where  b  -  2^— 
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By  measuring  the  amplitude  (9)  and  phase  angle  ,  (<p^)  of  the  angular 
velocity  9  at  different  frequencies  to  C  and  C  can  be  deter- 

in 

q  a 

mined  by  a  circle  or  Nyquist  plot  (see  reference  2),  where 
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tan 


-1 


U  ma  ,  I/P  SL3 
“L  C  U_/_  _  \ 

q.  tOL\  m  / 


S  (6) 


<Pe  (to)  +  | 


Therefore,  in  principle,  the  aerodynamic  derivatives  Cm  and 
C  can  be  found.  However,  it  will  be  shown  in  Section  2.1  that  this 

mq 

method  is  not  feasible  in  hypersonic  flow,  and  it  is  necessary  to  turn 
elsewhere  for  an  acceptable  technique. 


1.3  Free  Oscillation  Technique 

In  this  method  the  model  is  free  to  oscillate  in  the  sting 
bearings  as  before,  however  the  bearings  remain  stationary.  The  motion 
of  the  model  is  purely  oscillation  about  a  fixed  axis,  under  the  influence 
of  bearing  friction  and  aerodynamic  forces.  Here 


0 


w 

U 


a. 


w 

U 


a 


(7) 


4. 


and  Equation  (2)  becomes 

a  +  ba  +  ka  =  0  (8) 

Adding  an  arbitrary  bearing  friction  torque  f(6;)  gives  the  equation  of 
motion  of  the  system 

a  +  ba  +  f(a)  +  ka  »  o  (9) 


With  small  damping,  the  frequency  of  oscillation  is 
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C  can  be  simply  found  by  a  measurement  of  the  period  of  oscillation  of 

ma 

the  model,  and  determination  of  the  model  parameters  I,  S  and  L  and 


the  dynamic  pressure 


PU2 
2  ‘ 


The  derivative  C  is  a  measure  of  the  aerodynamic  damping 

mq 

present.  This  damping  in  turn  can  be  found  by  observing  the  rate  of 
decay  of  the  oscillation  amplitude  as  the  model  is  given  an  initial 
displacement  and  allowed  to  oscillate  in  the  flow  stream. 


2.  PRELIMINARY  CONSIDERATIONS 

From  the  beginning  of  this  program  it  was  apparent  that  the 
aerodynamic  damping  would  be  small  compared  to  the  restoring  moment. 
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and  using  estimated  values 

C  -  -  4/3  L  =  0.625  in 

°a 

c  -  -  2  U  =  2370  ft/sec 

°q 

then  b/k  ■  10~^/3,  an  extremely  small  amount  of  damping.  The  factor 
L/U  shows  the  difficulty;  in  hypersonic  flow  and  with  the  model  size 
restriction  imposed  by  the  flow  nozzle,'  L/u  becomes  very  small. 


2.1  Induced  Oscillation  Technique 
Equation  (6)  can  be  written 
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where  to  **  natural  or  undamped 

°  oscillatory  frequency 

-  VT 

Since  b  is  small,  co  will  be  very  close  to  and  any  slight 

fluctuation  in  this  difference  will  have  a  stroftg  effect  upon  the 
vector  position,  >/ 

A  preliminary  calculation,  using  topical  M  ™  6  low  density 
flow  values  and  a  45°  semi-vertex  angle  conpMtnOd’el  with  L  ■  0.625  in. 


6 . 


showed  that  a  frequency  co  =  ■ 1.0015  o>  would  shift  the  phase  vector  to 

cp  «  45°,  and  co  =  1.0060  <o  corresponds  to  cp  =  84.3°.  Tunnel  flow 

2 

characteristics  are  such  that  PU  fluctuates  as  much  as  ±  2%,  and  since 
coq  -VST.  coq  would  vary  ±  1%.  Such  a  natural  frequency  variation 
would  cause  almost  a  180°  oscillation  in  the  phase  vector,  rendering  it 
useless  in  determining  the  aerodynamic  coefficients  and  derivatives. 

Essentially,  the  small  amount  of  damping  forces  the  induced 

\ 

frequency  (co)  to  be  very  close  to  the  natural  frequency  (coq)  in 
order  to  use  the  Nyquist  plot,  and  with  the  flow  irregularities 
inherent  in  the  wind  tunnel,  coq  fluctuates  enough  to  make 
cp  =  tan  *  (coq2  -  co^)/cob  completely  undefinable.  This  conclusion 
resulted  in  abandonment  of  the  induced  oscillation  technique. 


2.2  Free  Oscillation  Technique 

Determination  of  C  by  observing  the  period  of  oscillation 

presented  no  serious  problems,  however  the  low  aerodynamic  damping  again 

gave  difficulty  in  finding  C  .  The  amplitude  decay  caused  by  bearing 

m 

q 

friction  would  be  relatively  large  and  could  easily  obscure  the  aero¬ 
dynamic  damping.  Information  on  low-friction  bearings  gave  a  minimum 
of  7  dyne-cm  of  resistance  torque  for  the  proposed  system.  Using  this 
figure,  calculations  at  M  =  6  flow  conditions  showed  that  the  damping 
from  the  bearings  would  be  approximately  10  times  the  aerodynamic  damping. 

On  the  basis  of  this  information,  it  was  decided  to  use  the 
free  oscillation  technique  to  determine  Cm  ,  with  the  possibility  that 

C  could  also  be  found.  C  would  be  a  function  of  the  difference  in 
m  m 

q  q 

decay  rate  between  the  model  oscillating  in  the  airflow  and  oscillating 
in  the  evacuated  tunnel  without  flow.  This  relationship  is  described  in 
detail  in  Appendix  A. 
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3.  EXPERIMENTAL  EQUIPMENT 

The  tests  were  conducted  in  the  No.  4  low  density  wind  tunnel 
at  the  University  of  California  Aeronautical  Sciences  Laboratory,  Richmond 
Field  Station.  The  Mach  6  nozzle  was  used,  providing  a  test  section 
approximately  2.5  inches  in  diameter.  The  fixture  holding  the  model  was 
mounted  directly  on  the  end  of  the  nozzle,  as  shown  in  Figure  1. 

Figure  2  shows  the  9°  cone  model  mounted  on  the  fixture.  The 
model  is  fixed  on  the  vertical  sting,  the  sting  being  supported  by  two 
electrically  vibrated  bearings.  A  cross  shaft  is  mounted  on  the  sting 
above  the  upper  bearing.  One  arm  of  this  shaft  carries  a  blade  which 
interrupts  the  light  beam  to  a  photocell  as  the  model  passes  through 
zero  angle  of  attack.  The  other  arm  is  attracted  by  an  electromagnet 
mounted  on  the  tunnel,  thereby  pulling  the  model  to  an  angle  of  attack. 
Release  of  the  arm  starts  the  oscillation.  In  addition,  a  shadow  cast 
by  this  arm  is  projected  upon  a  calibrated  screen,  providing  direct 
oscillation  angle  (ptQ)  information. 

As  shown  in  Figure  1,  a  cardboard  and  plastic  cover  shields 
the  arm  from  extraneous  air  currents.  The  quartz  fiber  is  shown  sus¬ 
pended  from  above  and  attaches  to  the  tip  of  the  sting.  The  sting  itself 
(here  carrying  the  45°  cone)  has  a  tungsten  carbide  center  section  2 
inches  long  and  0.078  inches  in  diameter.  This  design  represents  a  com¬ 
promise  between  sting  structural  stiffness  and  minimum  aerodynamic 
interference,  several  configurations  were  tried  with  this  being  the  most 
satisfac tory. 

The  photocell  output  is  fed  to  an  electronic  counter  which 
times  the  interval  between  pulses  and  presents  this  time  on  a  visual 
readout.  Since  the  swinging  arm  interrupts  the  photocell  beam  twice 
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complete  oscillation,  the  time  shown  is  one  half  the  period  of  oscilla¬ 
tion. 

The  bearings  have  an  electrically  vibrated  inner  bearing  race, 
which  reduces  friction  torque  to  a  minimum  by  eliminating  static  friction. 
They  are  the  most  suitable  bearings  that  could  be  found  for  this  applica¬ 
tion,  although  devices  such  as  air  and  magnetic  suspension  may  have  even 
less  friction.  These  possibilities  were  not  explored. 


4.  EXPERIMENTAL  PROCEDURE  AND  DATA  REDUCTION 


4.1 


The  model  was  given  an  initial  displacement  of  approximately 
35°  and  allowed  to  oscillate  in  the  flow  stream.  Oscillation  period  data 
from  the  electronic  counter  was  recorded  at  ClQ  =  30°,  20°,  15°,  10°, 

5°,  and  3°.  This  period  gave  C  (a  )  directly  through  the  relationship 


2 

21 

pu2Sl 


4.2 


A  stopwatch  was  used  to  determine  the  oscillation  decay  rate. 

As  the  motion  decayed  through  values  of  Q!0  ■  30°,  20®,  15°,  10°,  and  5°, 
the  time  was  recorded.  The  average  decay  rate  was  computed  in  the 
intervals  between  these  values  by  dividing  the  particular  interval  (in 
degrees)  by  the  time  the  model  took  to  decay  through  this  interval.  This 
value  was  taken  as  the  decay  rate  at  the  angle  in  the  mid-point  of  the 
interval.  For  example,  in  the  interval  between  15°  and  10°,  the  decay 
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rate  would  be  the  interval  (5°)  divided  by  the  decay  time  (in  seconds) 
from  15°  to  10°.  This  value  is  them  the  decay  rate  existing  at  aQ  ■ 
12.5°. 

As  described  in  Appendix  A,  this  decay  test  is  performed  under 

both  no  flow  and  flow  conditions*  the  difference  in  rate  of  decay  giving 

the  value  of  C  ,  .  Using  this  method,  Cw  can  be  found  for  a  ■  25°, 
m.  •  m  ,  o’ 

q 

17.5°,  12.5°,  and  7.5°. 

4.3  Other  Parameters 

The  moment  of  inertia  of  the  system  (I)  was  computed  as 
described  in  Appendix  B.  S  and  L  were  accurately  measured  using 
ordinary  machine  shop  techniques.  Flow  information  (P,  M,  U,  and  Re/in) 
were  computed  from  compressible  flow  charts  using  flow  pressure  and 
temperature  data  recorded  with  each  run. 


5.  RESULTS 

Two  cone  models  were  tested,  one  having  a  semi-vertex  angle 

7  **  45°,  and  the  other  7  ■  9°.  For  the  45°  cone,  x/l  ■  0.337,  for 

the  9°  cone,  x/L  ■  0.350.  The  average  flow  conditions  in  these  tests 

were  M  ■  6.05,  Re/in  ■  11,700,  U  ■  2375  ft/sec,  «. 0.0530  p$ii,  and  a 

static  pressure  of  125  microns  of  mercury. 

Values  of  C  are  tabulated  in  Figure  3  and  shown  in  graphi¬ 
te 

cal  form  in  Figure  4.  C ^  is  plotted  in  Figure  7,  These  values  are 

given  as  functions  of  o^,  the  angle  of  oscillation.  No  attempt  has 

been  made  to  define  them  as  functions  of  the  Instantaneous  angle  of 

attack,  (a),  and  since  both  C_  and  ,C^  vary  with  angle  of  attack, 

a  “q 

the  figures  given  represent  an  average  or  effective  value  over  a  single 
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oscillation.  Hereafter  these  effective  values  will  be  designated  as 

C'  and*  C.  .  . 
u  q 

6.  EXPERIMENTAL  ACCURACY 

A  numerical  error  analysis  is  givcrn  in  Appendix  C;  the  figures 

below  are  derived  therein. 

6.1  C., 

. .  ma 

The  uncertainty  in  the  values  of  C_  is  between  3.0%  and 

a 

3.6%,  depending  upon  aQ.  The  prime  contributors  to  this  inaccuracy  are 

U, .  T,  and  a  .  The  error  in  U  is  inherent  in  the  flow  stream  and 
o 

cannot  be  improved  with  the  present  tunnel  and  nozzle.  The  uncertainty 

in  T  is  based  upon  data  scatter;  refinement  of  technique  could  possibly 

lower  this  somewhat.  The  inaccuracy  in  CCq  is  discussed  fully  in  the 

next  section.  These  values  of  C  are  about  as  precise  as  presently 

ma 

possible;  it  would  be  difficult  to  better  them  significantly  without 
improving  U. 


The  small  magnitude  of  damping  and  the  relatively  crude  decay 
rate  measurements  lead  to  the  large  uncertainties  in  Cm>  .  It  is  felt 

3. 

that  a  more  sophisticated  experimental  system  is  necessai^  to  derive 
accurate  data. 

The  difficulties  arise  with  the  measurement  of  0£  and  with 

o 

the  bearing  friction.  The  present  method  of  reading  CtQ  (shadow  of 

swinging  arm  projected  on  a  calibrated  screen)  causes  the  29%  uncertainty 

in  decay  rate  data;  an  Improved  system  here  would  definitely  improve  the 

measured  value  of  C_  . 

m 

'  q 
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Aerodynamic  damping  is  small  compared  to  damping  from  bearing 
friction,  consequently  in  the  subtraction  of  decay  rate$  (see  Appendix  A) 
two  nearly  equal  numbers  with  moderate  uncertainties  etre  subtracted,  the 
difference  being  small  in  magnitude  and  having  a  large  uncertainty. 

A  solution  here  would  be  to  reduce  bearing  friction.  The 
bearings  used  are  claimed  to  be  the  best  conventional  bearings  available, 
however  devices  such  as  air  suspension  and  suspension  by  fibers  would 
warrant  investigation  for  further  experiments^ 

7.  DISCUSSION  OF  RESULTS 

7 . 1  Cm 
_ a 

Figure  4  presents  the  experimentally  derived  values  and  compares 
them  with  other  data,  both  theoretical  and  experimental. 

On  the  ordinate  a  »  0  the  small  angle  values  of  C  from 
o  m 

a 

3 

the  theory  of  Coakley,  Laitone,  and  Maas  and  also  from  Tobak  and 
4 

Wehrend  are  plotted.  They  are  seen  to  be  inreasonable  agreement  with 
the  experimental  values,  the  45°  cone  particularly.  Both  theoretical 
analyses  require  a  <  y  (all  elements  of  the  cone  must  "see"  or  face 
the  flow),  the  9°  cone  is  near  this  limit  even  at  small  0.q,  which  might 
explain  the  great  deviation. 

The  bottom  curve  is  from  the  experimental  data  of  Maslach  and 
Talbot^  (obtained  under  the  same  flow  conditions  as  this  report)  on  the 
static  moment  of  a  9°  cone.  From  their  plot  of  Cm  vs  a  (corrected 
for  center  of  oscillation  position)  the  slope  of  the  secant  line  between 
a  m  0  and  a  *  a  was  taken  as  the  average  slope  of  Cm  during  an 
oscillation  bounded  by  0tQ.  Therefore  the  values  plotted  are  Cm 


and 
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are  directly  comparable  to  the  experimental  values.  It  is  seen  that  C 

Of 

from  the  static  test  is  generally  smaller  than  that  from  the  dynamic  test 

and  varies  less  with  angle  of  attack.  These  differences  are  probably 

caused  by  unsteady  flow  effects  around  the  oscillating  cone  and  are  seen 

to  be  appreciable  in  this  hypersonic,  low  density  flow  regime. 

Charwat^  develops  a  Newtonian  flow  analysis  for  the  normal 

force  derivative  (which  i6  proportional  to  assuming  a 

stationary  center  of  pressure)  for  a  <  7  and  also  for  the  case  where 

half  of  the  surface  of  the  body  sees  the  flow.  For  a  cone  this  latter 

solution  is  strictly  valid  for  a  =  90°  only,  but  the  extension  of  this 

theory  to  lower  angles  of  attack  can  give  an  indication  of  the  variance 

of  C  with  a  (or  a  )  for  the  9°  cone. 

0 

The  two  solutions  of  Charwat  are  shown  in  Figure  5,  with  a 
curve  (dashed  line)  faired  between  the  two  being  more  representative  of 
the  actual  case.  It  is  seen  that  CN  (and  therefore  C  )  generally 

increases  with  angle  of  attack  up  to  a  «*  35°.  This  correlates  well 
with  the  experimental  values  of  C  for  the  9°  cone  shown  in  Figure 

_  "k 

4,  where  C  is  seen  to  increase  with  a  . 

Coakley,  et.  al.,  also  give  a  hypersonic,  Newtonian  flow 
motion  analysis  not  restricted  to  small  angles  of  attack,  requiring  only 
that  a  <  7.  Their  non-linear  equation  of  motion  is  (disregarding 
damping) : 

a  +  SL  ^  4/3  -  2  —  cos2  7  — -  0  (13) 

1  2  3 

Using  the  expansion  j  sin  2a  -  a  -  j  OC  +  ...  .  the  frequency  of 
oscillation  of  motions  governed  by  this  equation  is 


13. 


a)  -  2n  "\Jb£  ^i[4/3  -  2  |  cos2  7]  “0  -  -  ao2  (14) 

Figure  6  compares  this  relationship  p,  V,S,L  and  1  taken  the  same  as  in 

experiment)  with  the  frequencies  obtained  in  the  45°  cone  experiment. 

It  is  seen  that  the  decrease  in  03  with  larget  CtQ  is  even  stronger 

than  predicted  by  the  non-linear  theory. 

Apparently^  these  differences  are  caused  by  viscous  effects. 

Newtonian  flow  theory,  of  course,  does  not  consider  the  boundary  layer-. 

and  wake  phenomena  which  certainly  influence  the  motion,  especially  at 

large  a.  Very  little  is  known  about  hypersonic  viscosity  effects  upon 

dynamic  stability;  the  experimental  results  here  indicate  that  they  are 

appreciable  and  should  be  further  investigated^ 

It  should  be  noted  that  all  the  values  of  Cm  are  negative, 

a 

indicating  positive  static  stability  in  equations  (8)  and  (9).  For  the 

case  of  a  missile  returning  to  earth  through  the  atmosphere,  Allen^ 

shows  that  it  is  only  necessary  that  C  <0  for  static  stability, 

“a 

after  considering  the  effects  of  the  earth’s  gravitational  field, 
deceleration  due  to  aerodynamic  drag  and  variance  of  atmospheric 
density  with  the  altitude.  Thus,  it  seems  that  static  stability  of  a 
conical  missile  is  assured  as  long  as  the  center  of  gravity  (and  there¬ 
fore  the  axis  of  oscillation)- is  well  forward,  as  in  these  tests.  As 
Allen  points  out,  however,  static  instability  can  occur  with  a  c.g. 
position  towards  the  rear  of  the  missile. 

It  is  felt  that  the  values  of  C  found  in  these  experiments 
are  representative  of  the  models  tested.  It  is  not  possible  to  improve 
the  experimental  accuracy  significantly  without  improving  the  flow 


I 


14. 


if; 

I' 


characteristics  of  the  wind  tunnel.  However,  it  does  seem  that  tests 

with  the  center  of  oscillation  moved  aft  would  be  worthwhile,  since  a 

point  would  be  reached  where  C  >0  and  the  cone  is  statically  un- 

ma 

stable.  For  small  a  this  point  would  be  near  the  center  of  pressure, 
but  for  larger  oscillations  (as  in  these  tests)  the  point  of  instability 
cannot  be  predicted  theoretically,  and  a  wind  tunnel  test  is  called  for. 


7.2  C 


m 


The  values  of  C  shown  in  Figure  7  have  such  a  high  experi- 
mq 

mental  uncertainty  that  they  must  be  considered  in  a  qualitative  6ense 

only.  Very  little  can  be  said  of  their  absolute  magnitude,  except  that 

perhaps  C  for  the  A5°  cone  decreases  with  Of  and  for  the  9°  cone 
m  o 

q 

increases  with  a  . 

o 

It  should  be  realized  that  Cm  has  been  found  to  be  positive, 

*1  . 

indicating  negative  damping  and  dynamic  instability.  In  every  test  on 
both  cones,  the  model  took  longer  to  decay  (from  30°  to  5°)  when  under 
flow  than  it  did  with  no  flow.  This  is  indicative  of  energy  input  into 
the  model  from  the  airflow. 

This  instability  here  only  partially  counteracts  the  bearing 
friction  damping,  and  the  oscillations  still  decrease.  However,  in  an 
object  in  free  flight,  instability  would  tend  to  cause  increasing 
oscillations  (disregarding  the  influence  of  other  damping  factors,  of 
course).  For  this  reason  it  is  apparent  that  the  instability  found  here 
is  Important  and  should  be  investigated  further. 

All  the  available  theory  predicts  C_  <  0,  and  previous  ex- 
periments  (at  lower  Mach  numbers)  show  Qm  <  0  and  nearly  independent 

q 
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of  the  angle  of  attack.  It  must  be  emphasized,  though,  that  hypersonic 
Newtonian  flow  theory  does  not  consider  any  variation  in  the  pressures 

on  the  rear  of  a  cone,  and  this  effect  could  be  important  for  a  cone 
oscillating  in  a  flow  stream.  An  analytical  study,  one  including  viscous 
effects  and  not  restricted  to  small  a,  would  certainly  be  difficult, 
but  might  reveal  the  nature  of  this  instability. 

There  seem  to  be  no  prior  dynamic  stability  experiments  of 
this  type  in  hypersonic  flow.  Further  investigation  is  certainly  in 
order,  both  to  more  accurately  determine  C  ,  and  especially  to  investi- 

mq 

gate  the  nature  of  the  instability. 

As  stated  earlier,  the  large  uncertainty  in  Cm  stems  from 

mq 

the  small  amount  of  damping  present.  With  the  model  size  restriction 
of  the  f low  nozzle,  the  ratio  b/k  will  be  small  in  hypersonic  flow 
unless  C  (and  therefore  k)  in  reduced  by  oscillating  the  model  on 

mq 

an  axis  near  the  center  of  pressure.  Unfortunately,  reducing  the  restor¬ 
ing  moment  in  this  manner  increases  the  adverse  effect  of  bearing  friction, 
so  thi6  method  most  likely  would  not  be  an  improvement. 

Perhaps  a  free  flight  'shot'  in  a  ballistic-type  tunnel  would 
be  the  solution- -here  there  would  be  no  bearing  friction  and  the  model 
would  oscillate  about  its  center  of  mass.  High  speed  photographs  could 
provide  an  oscillation-time  history;  the  decay  (or  increase)  in  the 
oscillation  would  give  the  damping.  However,  deceleration  from  aerody¬ 
namic  drag  introduces  other  terms  into  the  damping  coefficient;  these 
might  prove  difficult  to  separate  from  Cm  . 

q_ 

More  accurate  determination  of  C  would  be  possible  by 

m 

q 

improvement  of  the  present  apparatus.  The  uncertainty  in  aQ  (±  1°)  is 
large  and  certainly  could  be  decreased.  An  improved  optical  system  or 
high  speed  motion  pictures  might  be  the  solution. 
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Decreasing  the  bearing  friction  would  increase  the  relative 
effect  of  the  aerodynamic  damping  and  thereby  improve  the  accuracy  of 
Cm  .  The  bearings  used  here  are  probably  as  friction  free  as  possible 

q 

in  a  mechanical  system,  but  perhaps  could  be  improved  slightly.  Magnetic 
or  air  suspension  might  be  feasible,  although  it  would  be  difficult  to 
avoid  the  air  or  magnetic  field  affecting  the  motion  of  the  model. 

Suspension  by  a  quartz  fiber  has  definite  possibilities,  for 
quartz  has  very  little  internal  friction  and  a  quartz  torsional  pendulum 
decays  very  slowly.  A  model  made  from  a  heavy  material  such  as  lead 
could  be  hung  in  the  airflow  from  a  quartz  fiber,  perhaps  with  another 
fiber  holding  it  from  below. 

Reducing  the  friction  of  the  suspension  device  and  decreasing 

the  uncertainty  in  aQ  would  permit  an  accurate  plot  of  vs  time. 

The  slope  of  this  curve  could  easily  be  found  by  use  of  a  least  squares 

program  on  a  digital  computer.  This  would  eliminate  the  averaging 

method  used  in  this  report  and  provide  accurate  values  of  decay  rate  at 

any  a  . 
o 

Another  avenue  for  further  investigation  is  the  cause  of  the 
negative  damping.  As  previously  mentioned,  there  has  been  no  prior  indi¬ 
cation  of  dynamic  instability  from  either  theoretical  or  experimental 
studies,  although  little  work  has  been  done  in  either  field  on  oscillating 
bodies  in  hypersonic  flow. 

Possibly  the  instability  is  due  to  the  asymmetrical  shedding  of 
vortices  from  the  sharp  trailing  edge  of  the  cone.  Very  little  is  under¬ 
stood  about  vortex  formation  on  three-dimensional  oscillating  bodies,  but 
it  is  possible  to  make  a  few  statements  concerning  their  effect. 


17. 

f 

It  is  known  that  vortices  depart  in  a  periodic  manner  from 
stationary  bodies  of  revolution  at  zero  angled of  attack,  and  that  they 
affect  the  pressure  distribution  on  the  surface'  of  the  body.  It  seems 
reasonable,  then,  that  pressures  on  an  oscillating  body  are  also  affected 
by  vortex  formation.  If  these  pressure  fluctuations  * occur  in  a  manner 
so  as  to  transfer  energy  from  the  flow the  model  (i.e.,  the  additional 

i 

torque  from  the  pressure  fluctuation  acts  in  the  same  direction  as  the 

>  ■ 

rotational  velocity,  rather  than  being  a . torque  opposing  the  velocity 

as  in  ordinary  positive  damping)  there  is  .negative  damping.  This  might 

well  be  what  is  happening  here. 

The  vibration  of  electric  power  lines  in  a  crosswind  is  well- 

known;  this  is  an  example  of  aerodynamic  instability  from  vortex  formation 

8 

on  a  two-dimensional  body.  Goldstein  speaks  of  the  oscillatory  nature 
of  vortex  departure  behind  a  flat  circular  disc  and  also  of  vortex  in¬ 
stability  behind  a  two-dimensional  cylinder.  Unfortunately,  none  of  the 
information  available  specifically  pertains  to  an  oscillating,  three- 
dimensional  body. 

As  a  feasibility  study,  this  program  was  successful  first  in 

that  the  Induced  oscillation  technique  was  shown  to  be  impractical,  and 

also  that  C  could  be  easily  determined  by  the  free  oscillation  tech- 
_  a 

nique.  C  was  found  to  a  rather  poor  degree  of  uncertainty,  but 

q 

improvement  of  the  experimental  apparatus  could  produce  more  accurate 


results 
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APPENDIX  A 
DETEEM IMAT I OH  OF  C 

m 


The  aerodynamic  derivative  C,  is  a  function  of  the  decay 

mq 

rate  of  the  oscillations  in  the  flow.  However,  bearing  friction  also 
contributes  to  decay;  this  effect  must  also  be  evaluated.  Two  separate 
decay  tests  were  made,  the  first  in  the  evacuated  tunnel  with  no  flow 
but  with  a  quartz  fiber  serving  as  a  torsion  spring.  The  second  test 
was  the  sapae  as  above  (quartz  spring  was  retained)  except  that  the 
model  was  exposed  to  the  airflow. 

The  equations  of  motion  for  these  systems  are: 

No  Flow 

a  +  da  +  fa  <=  0  (A-l) 

Flow 

a  +  ba  +  da  +  ka  +  fa  =  0 

(A-2) 

■  a  +  (bfd)a  +  (k+f)a  ■=  0 


where:  d  =  effective  damping  coefficient 
of  bearings  and  quartz  fiber 

f  *»  restoring  moment  coefficient 
of  fiber 


6  «*  logarithmic  decrement  of  oscillations 


6  » 


<“o>„ 


and  8  ■  ^  B 


where  B  ■  Of  coefficient  in  general 


I 
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(“o'n  "  ^o^n+l  T 
6  -  - g-j - -  jl 

o'n 

2  «Vn  '  <°oW 


B  - 


(a  ) 

'  o'n 


(a  ) 

on 


F  <°0>n  F 


D_  =  b  +  d 


NF  (Of  )  NF 
'  on 


D„„  -  d 


where  the  subscripts  F  and  NF  refer  to 
the  flow  and  no  flow  cases  respectively 


Since 


(b  +  d)  -  d 


(a  ) 

o'n 


[df  '  dnf] 


also  b.EjlL-  [-c ] 


Cm  <V  3 


41 


PUSL2 


DF  ~  PNF 

a 


(A- 3) 


at  any  aQ  can  be  found  by  measuring  the  decay  rates 
under  flow  and  no  flow  conditions.  This  analysis  assumes  that  bearing 
friction  produces  linear  (viscous)  damping.  The  bearing  friction  is 
most  .likely  nearer  sliding  or  Coulomb  friction,  but  here  this  non¬ 
linearity  is  neglected  and  has  been  replaced  by  the  "equivalent" 
viscous  damping. 
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APPEND IK  B 

MOMENT  OF  INERTIA  DETERMINATION 


The  mass  moment  of  inertia  (I)  was  found  by  oscillating  the 
system  in  an  evacuated  tunnel  using  the  quartz  torsion  spring  and  observ¬ 
ing  the  period  of  oscillation  on  the  counter.  Two  tests  were  made,  one 
with  the  system  as  used  in  the  flow  tests,  and  a  second  using  an  addi¬ 
tional  cross  arm,  of  known  moment  of  inertia  I',  attached  to  the  sting. 
Since 

T  ~  Il/2 

then 

T  /  j 

with  arm  /  w/arm 

ip  '  '  **  p 

without  arm  V  w/o  arm 


or 


*w/o  arm 


=  I 


v, 


arm 


w/o  arm 

r 

w/arm 


I' 


(!'  +  I) 


A 

w/o  arm 

p 

w/arm 


hL 


o  arm 


o  arm 


w/arm 


l 


(B-l) 


w/arm 


This  indirect  method  permits  accurate  determination  of  I 
without  knowing  the  spring  constant  of  the  quartz  spring.  This  test  was 
performed  with  both  the  9°  and  45°  cones,  and  in  both  cases  I'  31. 
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APPENDIX  C 
ERROR  ANALYSIS 


m. 


a 


The  following  under tain ty  values  have  been  assigned: 


X 

m 

±  2.07. 

\  / 

* 

I 

m 

±  0.47. 

*  * 

S  ,  L 

m 

±  0.27. 

* 

U 

- 

±  1.07. 

* 

T 

- 

1.07. 

* 

a 

±  17. 

J  1  2  *  * 

where  \-r  PU  ,  I  ,  indicate 


uncertainty 


uncertainty  in 
in  I,  etc. 

A  least  count  analysis  would  indicate  an  uncertainty  of  ±  0.87.  in 
* 


f  PU2 

l 

±  27,  is 


|  ;  however,  past  experience  has  shown  that  the  figure  of 
more  realistic  due  to  non-uniform  flow  in  the  nozzle. 
Using  the  4-m-s  method  of  error  calculation. 


m 


a 


V 


..  *,2  .  *.  2  1 1  2* 
(2T  r  +  (I  )  +  If  Pu 


*  <S*>2  +  a*)2 


-  ±  2.9% 


This  is  the  error  in  C  as  observed.  Since  there  is  an 

°a 

uncertainty  in  the  value  of  a  ,  an  additional  uncertainty  in  C  (a  ) 

O  ffl  o 

a 

arises.  Examples  for  both  models  will  be  given. 
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45°  Cone 


Cm  changes  about  0.009  per  degree  of  aQ 


a  ®  30° 
o 

0.009 

0.765 


1.27. 


ns 

a 


92  +  1.22 


3.17. 


a  -  15° 
o 

0.009 

0.896 


1.07. 


£  -  V2-9* +  ^ 

a 


3.07. 


9°  Cone 


C  changes  0.008  per  degree  of  a 

a  0 


a  ®  20° 
o 


0.008 

0.761 


-  1.17. 


m 


a 


2  2 

r  +  i.r 


3.17. 


aQ  ®  10°  (The  worst  case)  -  changes  0.014  per  degree  of  an 


0.014 

0.666 


-1*  °£a  -  V9- 


92  +  2.12  -  3.67. 


The  uncertainty  for  the  45°  cone  at  Of  «  12.5°  will  be  taken 

o 

for  Illustration. 

D  -  rat.  of  decay  of  ..dilation  -  BBStSU&a 

time  to  decay 

*qo>l  ~  (°o>2  15*  -  10c 

*2  "  fcl  t10  "  t15 
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The  scatter  of  data  indicates  that  the  average  time  figures  are  accurate 
within  ±  2%  ,  and  using  a*  =  ±1° 


is  a  direct  function  of  D_  -  D._  and  therefore  C,  will  have 

r  Nr  0 
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an  uncertainty  of  at  least  410%.  All  other  values  of  C 

ra 

( 

similarly  large  uncertainty,  this  example  being  typical. 


show  a 


FIG.  I  FIXTURE  MOUNTED  ON  NOZZLE  WITH  45°  CONE 

HYD  7610 


FIG.  2  FIXTURE  WITH  9°  CONE 
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9°  CONE 

3° 

1.007 

0.582 

5° 

0.994 

0.598 

10° 

0.949 

0.666 

1  5° 

0.896 

0.722 

o 

o 

CM 

0.851 

0.762 

o 
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ro 

0.765 

0.803 

FIG.  3  VALUES  OF  Cmo(  FOR  VARIOUS 
OSCILLATION  ANGLES  (<*0)  . 
45°  t,  9°  SEMI- VERTEX  ANGLE 
CONES. 


HYD  7612 


HYO  7613 
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